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ABSTRACT 

The  stability  of  a  boundary  layer  on  a  heated  flat  plate  is  investigated  in  the  linear 
regime.  The  flow  is  shown  to  be  unstable  to  longitudinal  vortex  structures  which  in  general 
develop  in  a  nonparallel  manner  in  the  streamwise  direction.  Solutions  of  the  nonparallel 
equations  are  obtained  numerically  at  0(1)  values  of  the  appropriate  stability  parameter, 
ie  the  Grashof  number.  The  particular  cases  investigated  relate  to  the  situations  when  the 
instability  is  induced  by  localized  or  distributed  wall  roughness  or  nonuniform  wall  heating. 
The  case  when  the  vortices  are  induced  by  freestream  disturbances  is  also  considered.  The 
fastest  growing  mode  is  found  to  be  governed  by  a  quasi-parallel  theory  at  high  wavenumbers. 
The  wavenumber  and  growth  rate  of  the  fastest  growing  mode  are  found  in  closed  form.  At 
low  wavenumbers  the  vortex  instability  is  shown  to  be  closely  related  to  Tollmein-Schlichting 
waves,  the  effect  of  wall  heating  or  cooling  on  the  latter  type  of  instability  is  discussed. 
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1  Introduction 


Our  concern  is  with  the  insta]:)ility  of  forced  convection  boundary  layers  over  hori¬ 
zontal  heated  flat  plates.  Such  flows  are  unstable  to  at  least  two  types  of  hydrodynamic 
instabilities;  firstly  we  expect  a  convective  Rayleigh- Benard  type  of  instability  because 
the  fluid  at  the  plate  is  hotter  than  the  fluid  in  the  freestream.  Secondly  we  expect  a 
Tollmien-Schlichting  type  of  instability  because  of  the  similarity  of  the  flow  to  isothermal 
i  boundary  layers  where  that  type  of  disturljance  is  known  to  be  important.  In  this  paper 

we  shall  be  primarily  be  concerned  with  the  vortex  mode  of  instability  which  we  iiives- 

I 
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tigate  in  a  self-consistent  manner  using  an  approach  suggested  ])y  related  w'ork  on  the 
closely  connected  Gbrtler  vortex  problem,  see  for  example  Hall  (1990)  and  Denier,  Hall 
and  Seddougui  (1991).  However  in  the  small  wavenumber  limit  of  the  vortex  mode  we 
find  an  unexpected  relationship  between  the  vortex  mode  and  Tollmien-Schlichting  waves; 
in  effect  we  find  that  in  this  limit  the  two  modes  coalesce.  We  are  thus  able  to  describe 
both  propagating  vortex  modes  and  determine  the  effect  of  wall  heating  on  the  growth  of 
longitudinal  vortex  structures  in  boundary  Layers. 

Interest  in  forced  convection  boundary  layers  is  generated  by  the  wide  range  of  prac¬ 
tical  problem  where  such  flows  occur;  in  pju  ticular  we  refer  to  the  heat  transfer  problems 
associated  with  solar  heating,  electronic  devices  and  nuclear  reactors.  In  such  situations 
it  is  important  to  know  the  parameter  regime  where  instability  begins  because  of  the 
associated  change  in  heat  transfer  properties  of  the  flow. 

Experimental  investigations  of  tlu*  vort('x  mode  of  instability  in  a  forced  convection 
boundary  layer  have  been  carried  out  1)y  Gilpin.  Imura  and  Cheng  (1978)  for  water,  and 
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by  Wang  (1982)  for  air.  Here  Ave  shall  concentrate  on  boundary  layers  in  air  and  attempt 
to  explain  some  of  the  results  found  by  Wang.  In  fact  both  authors  demonstrated  the 
existence  of  the  onset  of  a  vortex  mode  of  instability  and  suggested  that  this  onset  occurs 

—  3 

at  the  same  value  of  ^  when  the  local  Grashof  number,  G^,  and  Reynolds  number 

Rx  are  varied.  We  also  note  that  a  relat('d  instability  occurs  in  channel  flows  when  one 
wall  is  heated,  the  reader  is  referred  to  the  paper  by  Akiyama,  Hwang,  and  Cheng  (1971) 
for  an  experimental  investigation  of  that  problem. 

Theoretical  investigations  of  the  vortex  mode  of  instability  have  been  given  by  Wu 
and  Cheng  (1976)  and  Moutsoglou,  Chen  and  Cheng  (1984).  In  both  of  these  calculations, 
and  all  other  investigations  we  are  aware  of,  the  growth  of  the  boundary  layer  is  not  taken 
into  account  in  a  self-consistent  manner.  In  effect  it  has  been  previously  assumed  that 
the  streamwise  variation  of  the  vortex  mode  is  on  a  short  Icngthscale  compared  to  that 
over  which  the  basic  flow  evolves.  However,  at  finite  values  of  the  Grashof  number,  where 
instability  first  sets  in,  the  destabilizing  buoyancy  forces  are  sufficient  only  to  provoke  a 
response  on  the  same  streamwise  lengthscale  as  that  over  which  the  basic  state  develops.  In 
that  case  it  follows  that  previous  calculations  have  ignored  a  crucial  property  of  the  vortex 
instability;  thus  the  intrinsic  nonx^aralh'!  nature  of  the  disturbance  has  been  neglected. 

A  similar  parallel  flow  assumption  was  made  for  many  years  by  researchers  concerned 
with  the  Gortlcr  mechanism  in  curved  l)oundary  layers.  More  recently  Hall  (1983)  showed 
that  the  nonparallcl  nature  of  the  vortex  mode  at  finite  values  of  the  Gortler  number 
must  be  accounted  for  by  numerical  investigations  of  the  disturbance  equations.  At  high 
Gortler  numbers  analytical  iirogress  can  be  made  because  the  growth  rates  become  large 


and  nonparallel  effects  may  be  neglected  at  zeroth  order;  sec  Hall(  1982a, 1982b),  Denier 
Hall  and  Seddougni  (1991). 

In  this  paper  we  shall  concentrate  on  tlu'  vortex  mode  at  finite  values  of  the  appropriate 
stability  parameter,  ie  the  Grashof  numbon’.  A  que.stion  of  fundamental  importance  in 
this,  and  in  fact  any  other  parameter  regime,  is  that  of  what  is  the  physical  process  which 
triggers  disturbances  in  the  flow  and  causes  tla'iii  to  amplify.  This,  the  so-called  receptivity 
problem,  has  been  addressed  for  the  Gdrtler  problem  by  Hall  (1990),  and  Denier,  Hall  and 
Seddougui  (1991).  In  particular  the  receptivity  problems  for  wall  roughness  and  freestream 
disturbances  wei'c  investigated  in  those  papers.  The  corresponding  receptivity  problems 
will  be  investigated  here  for  heated  boundary  layers. 

Thus,  following  a  formulation  of  the  instability  equations  in  Section  2,  we  shall  in 
Section  3  investigate  the  stability  problem  at  finite  Grashof  numbers  using  a  numerical 
scheme  outlined  in  Section  3.  In  Section  4  we  c<niceutrate  on  the  amplification  of  vortex 
disturbances  induced  by  localized  wall  roughu<.“ss  or  nonuniform  wall  heating.  In  Section 
5  we  will  investigate  the  generation  of  v<jrtex  disturbances  by  freestream  inhomogeneities, 
whilst  in  Section  6  some  results  for  distril)utcd  forcing  are  given. 

In  Section  7  we  concentrate  on  the  high  Grasluif  number  limit,  this  regime  is  relevant 
to  disturbances  which  have  passed  through  the  order  one  Grashof  numljer  regime  without 
sufficient  amplification  to  lx;  coiitrolh'd  In-  nonlinear  effects.  At  order  one  wavenumbers 
we  show  that  the  now  unique  growth  rate  is  determined  by  an  inviscid  eigenvalue  problem. 
However  the  growth  rate  predicted  Ijy  Hie  in\-iscid  theory  increases  monotonically  with  the 
wavenumber  so  that  the  fastest  growing  mode  cannot  be  descrilied  by  a  purely  inviscid 
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theory.  In  fact  at  sufficiently  high  wavemunbers  viscous  effects  come  back  into  play  and 
the  fastest  growing  mode  is  determined  a  viscous  eigenvalue  problem.  At  small  vortex 
wavenumbers  viscous  effects  again  dominate  and  we  show  that  the  vortex  mode  ultimately 
takes  on  a  triple- deck  like  structure.  This  structure  is  shown  to  describe  both  vortex  and 
Tolhnien-Schlichting  modes.  Finally  in  Section  8  we  compare  our  I’esults  with  previous 
experimental  and  theoretical  work  and  draw  some  conclusions. 

2  Formulation  of  the  Instability  Equations 

We  consider  the  flow  of  a  viscous  fluid  ov<'r  a  heated  semi-infinite  flat  plate.  Suppose 
that  t/oo  is  a  typical  freestream  velocity.  T  is  a  ty2:)ical  lengthscale  in  the  streamwise 
direction  and  v  is  the  kinematic  viscosity.  Th<'  Reynolds  number  is  defined  by  Re  — 
and  throughtout  we  assume  that  Re  »  1.  The  wall  is  defined  with  resi^ect  to  dimensional 
Cartesian  coordinates  x*,y*,z*  by 

y*  =:  LARc-'^f{x*IL,Re-^z*IL) 

where  A  is  a  small  dimensionless  constant.  We  take  the  temperature  of  the  fluid  a  long 
way  from  the  plate  to  be  zero  whilst  at  the  ])late  the  tcmi^erature  is  given  by 

T*  =To{l+A<j(,-*/L.Rc'^z*/L)). 

Here  Tq  is  a  constant  reference  temperature  and  g  represents  the  effect  of  a  slight  nonuni¬ 
form  heating  of  the  plate. 

We  define  (x,ij,z)  by 

ix,y.r.)  =  (x*/L.nAy*IL.nc'^z*IL) 
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and  a  dimensionless  velocity  vector  by  writing 

,  4-  -I-  +.  {u*,Re^v*,Re^w*) 


(2.1a) 


We  take  the  corresponding  pressure  function  to  be 


p+  =  p(x)  +  ARe  ^p(.i-,y,z)  +  0{A^)  = 

^  ooP 


(2.1t) 


where  p  is  the  fluid  density  and  p*  is  the  dimensional  fluid  pressure.  We  now  write 

(u'*’,  i;'*’,  w'^)  =  (77,  r,0)  +  A(i7,  C\72:)  +  O(A^),  (2.1c) 


and  the  temperature  field  then  expands  as 
T*  , 


=  T+  =  Re--^T{x,y)  +  ARe-H{x,y.,z)  +  OiA'^)  . 


(2.1d) 


In  the  forced  problem  for  the  order  A  field  the  basic  velocity  and  temperature  fields  u,  v,  T 
will  be  known  functions  of  x  and  y.  The  perturbed  velocity  and  temperature  fields  u,  v,  w,  0 
depend  on  all  three  dimensionless  coordinates.  The  steady  Navier  Stokes  equations  for  the 
problem  involving  buoyancy  forces  are 


u*u*.  +  v*u*. 


u*v*.  +  +  tu*vZ.  - 


1  dp* 
P  dx* 


1  dp* 


(2.2a) 

:•  ‘-’y*  y-  +  z*  )> 

(2.26), 

(2.2c) 

where  g  is  the  acceleration  due  to  gravity  and  is  the  coefficient  of  expansion.  The  energy 
equation  takes  th*'  form 


n*T*  +  v*T*  +  w*T:, 

^  y  - 


(t;.,. pt; 

pCp 


(2.2c/) 


o 


where  K  is  the  thermal  conductivity  and  Cp  is  the  constant  specific  heat.  It  is  to  be  noted 
that,  since  we  are  considering  a  low  speed  subsonic  flow,  the  term  representing  the  net  rate 
at  which  shearing  forces  perform  work  on  the  flviicl  is  negligible. 

We  now  write  (2.2a-d)  in  dimensionless  form  and  consider  the  limiting  form  of  these 
ecpiations  when  A  is  small  and  the  Reynolds  number  is  large. 

At  leading  order  order  we  obtain 


a  X  ^  y  —  0 , 

(2.3«) 

Tiu-x  "b  y  —  Px  “b 

(2.36) 

1 

UTx  +  VTy  =  -Tyy, 

a 

(2.3c) 

for  the  basic  velocity  and  temperature  fields  where  a  is  the  Prandtl  number.  At  next  order 
we  obtain 


Vx  +  I’v  +  =  0,  (2.4a) 

uuj:  4-  uiix  +  vHy  +  vhy  =  Uyy  +  u..,  (2.46) 

UVx  +  UVx  +  VVy  +  i'Vy  =  -fjy  +  Vyy  +  U.j  +  Gd,  (2.4c) 

iiiOx  +  ed'j,  =  -p:  +  Wyy  +  u>j-,  (2  Ad) 

u9x  +  ^9y  +  ufx  +  vfy  =  -{9yy  +  Ozt)-,  (2.4e) 

(7 


where  G  is  the  Grashof  number  defined  l>y 

G  =  L'oTodR^ 


G 


The  basic  flow  is  taken  to  be  the  Blasius  velocity  profile 

«  =  /'('/).  =  -^(?//'  -  /), 

V2x 

where  /"'  +  //"  =  0  with  /(O)  =  /'(O)  ==  0, 

/'(0G)  =  1. 

Here  the  similarity  variable  rj  is  defined  b}' 


The  basic  pressure  p  is  then  zero  and  the  liasic  temperature  profile  is  then  a  function  of 
just  77.  The  boundary  conditions  on  the  order  A  field  are  obtained  by  setting  up  Taylor 
series  expansions  about  7/  =  0;  we  obtain 


U  —  -fUy,V  = 

=  w  ■■ 

=  0,  0  =  g  -  ffy,  y  =  0,  ^ 

(2.5a) 

u,  V,  w,  9  0, 

V 

—*•00,  j 

U  =  f)  = 

0, 

r  =  l,  y  =  0  ) 

) . 
t . 

(2.56) 

u  — 7  1, 

T  - 

0,  y  — >•  00  J 

We  assume  that  /, </  are  such  that  we  can  write  fiiy  =  —Fi{x)q{z),  g  —  fTy  = 
F2{x)Q{z)  in  (2.5)  which  enables  us  to  Fourier  transform  the  disturbance  equations  in  the 
2  variable.  Thus  if  we  Foirrier  transform  the  cqufitions  (2. 4), (2. 5)  with  a  as  the  transform 
variable,  and  denote  the  transform  of  r],  1'/.  etc  Iw  q.ii  etc.  we  obtain 


i't  k  <’//  +  ><>«'  ~  0,  (2.6o) 

UUj:  +  VUy  +  VUy  +  U(/j  =  Uyy  —  U  ,  (2.66) 

uv^  +  vvy  +  vVy  +  nvr  —  -Py  T  Vyy  -  (1^ +  GO,  (2.6c) 

Xl.Wr  +  X’Wy  =  —p,  +  Wyy  -f  IC;;,  {2.6(1) 

ii,Tj:  +  uOj.  +  rTy  +  X'Oy  =  ~{0yy  —  (i~0),  (2.6f ) 
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Hence,  given  the  basic  velocity  and  temperatiue  profiles,  then  by  solving  the  system  of 
equations  (2.G6),  (2.6e),  (2.7)  subject  to  the  boundary  conditions  (2.C/)  we  can  determine 
the  solution  of  the  forced  convection  problem  by  numerical  integration  for  -finite  values  of 
the  Grasliof  number.  In  the  following  section  we  shall  describe  a  numerical  scheme  which 
we  have  used  to  integrate  the  disturbance  equations  found  above. 


3  The  numerical  scheme  and  some  preliminary  results 

The  disturbance  equations  (2.6b),  (2.Ge),  (2.7)  are  parabolic  in  x  and  hence,  having 
imposed  an  initial  disturbance  on  the  flow,  we  can  march  the  equations  downstream  from 
the  position  where  the  forcing  Ijcgins  and  monitor  the  vortex  growth  or  decay.  The  par¬ 
tial  differential  equations  describing  the  p<n'turbed  velocity  and  temperature  fields  were 
integrated  using  a  spectral  collocation  method  with  riliebychev  polynomials  used  to  ap¬ 
proximate  the  normal  dependence  of  the  disturljance.  The  Chebychev  polynomials  are 
defined  on  Q  e[— 1, 1]  by 


Tr(o)  =  cos(/.'cos  '  o  )  for  k  =  0,1,2 .. . 


S 


and  wo  approximate  the  streamwise  velocity  component  u,  for  example,  using  an  nth  degree 
Chcbychev  })olynoimal 


„(o)  =  y  =  2,7,.  =  =  1,2,... 

k=o 

and  are  the  unknown  Chcbychev  coefficic'iits  of  the  expansion.  The  change  of  variable 

(x,  y)  to  (a-,  y)  where  y  =  is  made  in  the  disturbance  ecpiations  which  are  then  solved 

\/2x 

on  r?e[0,  ?;oo]-  The  n  +  1  collocation  points  were  chosen  to  be 


a;  =  cos 


[i  -  1)t 


,  1  <  i  <  n  +  1. 


We  denote  the  value  of  u  at  y  =  ?/,,  xj  =  i'  +  je  by  uj  where  yi  =  -^(o/i  +  1),  x 
is  the  position  at  which  the  initial  disturbance  is  imposed,  J  is  an  integer  value  and  e 
is  the  steplength  in  the  x  direction.  A  similar  notation  was  employed  for  the  other  flow 
quantities.  We  suppose  these  cpiantities  are  known  at  the  Jth  step  and  illustrate  how  they 
are  advanced  to  a:  =  x^j+i  =  x  +  (j  +  l)e.  Consider  the  .r-momentum  equation  (2.Gb): 


J  /o  ^  -3  -3 


where  has  been  replaced  by  it'^  tinit<'  differc'iicc  ai)proximation 


ft?"*" '  —  ui 


We  approximate  by 


y  ^r,(o), 
t  =:0 
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and 


are  able  to  generate  the  successive  Cliebycliev  polynomials  using  the  relationship 

T, =  2zn(  =  )  ~  rt_,(c)  for  k-  >  0. 

In  a  similar  manner  we  can  determine  a  relationship  for  the  dmivatives  of  the  Chebychev 
polynomials.  We  rewrite  (3.1)  in  the  form 


for  1  <  —  1  with  f  =  1,  n  -f  1  for  each  value  of  /,'. 

We  use  the  streamwise  velocity  boundary  conditions  (2.G/)  to  replace  the  k  —  0,  k  =  n 
values  in  (3.2).  The  righthand  side  of  (3.2)  is  known  and  by  inverting  the  square  full  matrix 
on  the  left  hand  side  using  a  Gauss-.Jordaneliniination  inethocl,  the  values  of  the  Chebychev 
coefficients  for  0  <  A-  <  n  can  be  detenuiued.  In  a  similar  manner  we  can  update  (2.7) 
and  (2.Ge).  The  method  is  totally  implicit  so  that  we  c-xpect  to  have  a  numerically  stal^le 
scheme  for  a  streamwise  steplength  comparabh'  with  the  vertical  steplength.  The  spanwise 
velocity  component  w  can  then  be  calculated  from  the  contimiity  equation.  The  parameter 
7oo  was  varied  and  i/r^  =  10  was  found  to  provide  sufficiently  accurate  results.  The  number 
of  collocation  points  n  in  the  r/  direction  was  clK>sen  to  be  ^20  and  the  Prandtl  number 
(7  was  chosen  to  be  0.72.  the  value  for  air.  The  calculations  were  carried  out  on  an  .AMT 
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DAP510  and  the  code  was  written  so  as  to  take  advantage  of  the  architecture  of  that 
machine. 

In  order  to  validate  our  scheme  we  carried  out  some  calculations  for  the  case  when  the 
initial  form  of  the  disturbance  is  inipo.s('d  and  not  provided  by  a  receptivity  calculation. 
In  that  case  the  boundary  conditions  at  the  plate  and  far  from  the  wall  are 


u:=v~  —  =  6-  0,  y  =  0, 

oy 


=  —  =  0  =  0,  y  =  oo, 

oy 


(3.3a) 

(3.36) 


We  also  require 


n  =  UB(y)  ,  u  =  tiyly)  ,  0  =  SB(y),  x  =  x,  (3.4) 


where  the  initial  conditions  (3.4)  descril^e  some  given  vortex  perturbation  imposed  on  the 
flow  at  X  =:  X.  This  initial  disturbance  must  lie  consistent  with  (3.3).  Further  constraints 
on  the  initial  perturbation  (3.4  )  arc  required  in  order  to  avoid  singularities  in  the  velocity 
and  temperature  fields  at  t  =  J,  y  =  0 .  If  we  expand  u.t\9  as  Taylor  scries  about  x  =  'x 
and  y  =  0  we  find  that  the  required  conditions  are: 

a'/j (0)  =  0  ,  u';;(0)  =  a\d^(0),  (3.5a) 

u))’(0)  =  2a^/;)(0),  .(3.5b) 

e'l,{0)  =  0,  0"'(O)  =  a20'«(O).  (3.5c) 


The  perturbation  imposed  on  tin*  flow  was  taken  to  be  cither 


v/i  =  0  .  Oi)  =  0  ,  X  =  X, 


11 


or 


It  13  —  0  ,  VB  =  0  ,  (^B  ~  ^  ^ 


where  77  =  ——=  .  Note  that  both  (3.C)  and  (3  7)  satisfy  the  conditions  (3.3),  (3.5).  The 
\/2x 

numerical  scheme  described  in  this  section  was  used  to  solve  the  linear  disturbance  equa¬ 
tions  (2.6b),  (2.Gc),  (2.7),  subject  to  th  '  boundary  conditions  (3.3)  with  initial  conditions 
given  by  cither  (3.C)  or  (3.7). 

For  the  calculations  reported  here  we  took  e  =  0.1.  Of  course  the  accuracy  of  our 
calculations  was  checked  by  varying  e  and  the  vertical  grid  spacing  in  some  cases.  The 
vortex  growth  downstream  is  determined  by  monitoring 

£.  =  rwd'j 

Jo 

1  ^ 

and  the  local  growth  rate  (Ti  (x)  =  — —  .  The  neutral  point  was  taken  to  the  downstream 

Jhi  o.v 

location  where  this  growth  rate  vanishes,  the  local  Grashof  and  wavenumber  corresponding 
to  this  point  are  then  obtained  from 


=  a.r’/-  ,  G.  =  Gx^/^  . 


Different  neutral  curves  were  generated  for  fixed  G  =  0.025  by  varying  the  location  x 
of  the  initial  disturbance.  Figures  (3.1a-d)  demonstrate  the  downstream  velocity  and 
tem])eratnre  fields  for  an  initial  disturljance  given  by  (3.4)  with  G  =  0.025  and  a  =  0.069. 
The  corresponding  neutral  curves  are  given  in  Figure  (3.2).  For  an  initial  disturbance  of 
the  form  (3.7),  the  development  of  the  velocity  and  temperature  fields  is  shown  in  Figures 
(3.3a-d)  with  the  corrc'sponding  lu'utral  curves  shown  in  Figure  (3.4).  The  perturbed 
velocity  profiles  for  th('  initial  distnrl)auce  given  by  (3.4)  is  shown  in  Figures  (3.1a-c) 
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take  the  same  form  as  those  found  by  Hall  (1983).  The  essential  shape  of  the  perturbed 
velocity  and  temperature  components  do  not  alter  radically  with  increasing  x.  Initially  the 
spanwise  velocity  component  iv  is  proportional  to  du/dx  from  continuity,  since  the  normal 
velocity  component  v  of  the  initial  disturbance  is  zero.  Hence  both  the  streaimvise  velocity 
component  u  and  w  initially  decay  downstream  of  x  before  growing  at  larger  values  of  x. 
The  velocit}'^  and  temperature  profiles  for  a  perturbation  descriljed  by  (3.7)  are  of  a  similar 
form  but  take  the  opposite  sign.  It  is  clear  from  Figures  (3.2),  (3.4)  that  the  concept  of 
a  unique  neutral  curve  is  untenable  and  that  the  growth  or  decay  of  the  resulting  vortex 
structure  is  dependent  upon  its  initial  form  and  location.  However  far  downstream  the 
flow  is  disturbed  the  growth  rate  <7i(.r)  is  initially  negative.  On  the  right  hand  branch  of 
the  neutral  curve  Gi  ~  Uj  *,  but  for  a  fixed  wavenumber  disturbance,  Gx  ~  as  the  flow 
develops  >  0  and  the  flow  is  locally  stable.  Thus  any  disturbance  of  fixed  wavelength 

will  ultimately  be  stable  sufficiently  far  downstream  of  the  leading  edge. 

4  The  localized  forcing  problem 

We  now  consider  the  case  when  the  wall  forcing  is  described  by  isolated  forcing  func¬ 
tions  and  therefore  allow  Fi{x),  F^ix)  to  vary  on  a  relatively  fe  •  0(e)  lengthscale,  we 
write 

A'  (4.1) 

e 

where  e  is  small  and  we  assume  that  the  forcing  starts  at  x  =  x.  By  employing  a  faster 
streamwisc  lengthscale  for  the  isolated  heating  function  probh'in  we  will  provide  unique 
initial  conditions  for  the  disturbanc<- equations.  We  take  Fj,2(-'')  =  choose  x  = 

f  without  loss  of  generality.  This  fix<'s  the  original  lengthscah-  L  in  terms  of  the  distance 
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between  the  leading  edge  and  the  starting  pt)iiit  the  forcing.  In  the  following  discussion 
we  have  also  taken  q  —  Q  =  the  more  general  case  can  be  recovered  by  inserting  these 
factors  in  our  final  results.  In  order  to  find  the  forced  flow  in  a  neighbourhood  of  x  we  note 
that  uux  ~  Uyy  for  small  y  if  ?/  ~  0(e^/^)  and  hence  the  convection  and  diffusion  effects 
are  of  the  same  order  of  magnitude  in  a  layer  of  depth  For  small  y  the  basic  velocity 
and  temperature  fields  can  be  approximated  by 

{u,v,T)  =  (Ay,/ty-,l  +  ^ij)  +  ■  ■  ■ , 

where  A  =  u^(^,0),  —  Vyy{  0)  (p  =  Ty{  ~,0)  and  the  wall  forcing  implies  that  u,  6  are 
0(1)  in  the  region  y  ~  0(e^/^).  We  define  ^  and  assume  that  ^  =  0(1).  The 

appropriate  expansions  near  the  wall  are  then  found  to  be 

(u,  y,  u>,p,  d)  =  (no(A',0  +  •  •  •  0  +  •  •  • ,  €~'^Wo{X,  0  +  •  •  •  > 

e5/Vo(A-,0  +  .  ..>^0(^,0  +  ...)  .  (4.2) 
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where  Uo(5,^)  is  the  Laplace  transform  of  i)o(A",^),  Ai  is  the  Airy  function  and  A  and  B 
are  constants.  Transforming  the  continuity  equation  and  evaluating  it  at  the  wall  yields 


D  =  -sFx{s), 


where  F\  (5)  is  the  Laplace  transform  of  A’  )  .  Furthermore  we  require  that  Vo^  vanishes 


at  infinity,  so  that 


Vo  =  -sF(s)<  ^ 


Ai{y)dxj 


The  transformed  temperature  perturbation  and  streamwise  velocity  component  are  ob¬ 
tained  from  (2.6b,e)  which  give 


^  -  A^.s )  Uo  =  OoA, 


1^2  Ax 

--777  -  A^.s  jdo  =  VoX, 
a  d^2  y 


where  Qo  is  the  transformed  zeroth  order  temperature  field.  For  large  ^  the  asymptotic 
forms  of  U01V0  ai:id  9o  are  given  by 


\o  ~  +  .  ■  •  =  +  . . .  (4.Sa) 


)o^‘iFxs^l^\-^l^w  +  ... 


(4.86) 


So  ~  q-  . . . 


(4.8c) 
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where  u  =  — y4i'(0).  Hence  the  flow  within  th"  wall  layer  region  induces  the  motion  of  the 


fluid  in  the  y  —  0(1)  region  where  are  expanded  as 

»  =  (4.9a) 

a  =  +  (4.9fc) 

^  =  ^  (4  9c) 

and  the  Laplace  transforms  of  Uc  Vq  and  ©o  arc  found  to  be  given  by 

~  F  u' 

Uo  =  -3-4-5“‘/"A-’/Vm(y,a)  +  . . .  (4.10a) 

u 

Vo  =  3Fis2/3A-'/^u.’m(y,a)  +  . . .  (4.106) 

00  =  —3  l_  X~‘*^^u>m(y,  a)  +  . . .  (4.10c) 

u 

where  m  is  a  solution  of  the  stationary  Rayleigh  equation  problem. 

u{  dj^  —  a^)m  —  u"m  =  0,  (4.11a) 

7n(0)  =  1  ,  m(oo)  =  0.  (4.116). 


The  functions  in  (4.10)  decay  to  zero  cx2)onentially  as  y  — >  oo  and  satisfy  the  matching 
conditions  (4.8).  We  can  invert  the  Laplace  transforms  and  use  the  large  A'  form  of  the 
velocity  and  temperature  fields  as  the  initial  conditions  for  the  solution  of  the  perturbation 
equations  (2.Cb),  (2.6e),  (2.7)  subject  to  (2.Cf).  For  the  isolated  wall  functions  = 

6(A)  inversion  of  (4.5)  and  the  corresijonding  forms  for  Uo  and  9o  yields  the  similarity 
solution 


where  Ua,Va,Oa  are  known  functions  of  . 

A”  * ' 

We  now  demonstrate  how  the  similarity  solution  can  he  obtained  directly  from  the 
disturbance  equations.  The  similarity  variable  is  chosen  to  be  ^  —  y  j  —  |  where  x  = 


X  —  X.  For  (^  =  (9(1)  we  expand  the  pertuiiied  velocity  and  temperature  components  as 


rn  I  I 


(4.13a) 


(4.136) 


(4.13c) 


and  impose  the  wall  conditions 


uo  =  u;  =  uo  =  =  0 ,  ,e  =  0  . 


(4.13(i) 


The  above  expansions  are  substituted  into  (2.7)  and  at  leading  order  yield 


i,""  -I-  -I-  ==  0 


(4.14) 


where  a  prime  denotes  differentiation  with  re.spect  lo  if.  Equation  (4.14)  can  be  solved 
for  v'J  in  terms  of  Whittaker  functions  and  it  can  be  shown  that  there  exists  a  solution 
of  (4.14)  for  winch  Vg  =  v'^  =  9  ,  f  =  0,  and  such  that  at  infinity  ~  1+  exponentially 
small  terms.  It  is  to  be  noted  that  the  solutions  of  (4.14)  winch,  for  large  f,  behave  like 
f  and  terms  do  not  apiiear  in  the  required  solution  for  Vg.  Indeed  if  these  terms  did 
occur  the  wall  layi’r  solutions  could  not  lie  matched  in  the  rj  =  (9(1)  region  of  the  flow. 
Equations  (2.61)),  (2.Ge)  yield  at  leading  order 


II  ^  I  ,  c  _ 

>‘o  +  ~  ~\' 


(4.15a) 
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(4.156) 


-«o  +  3-«.  -f  («o  -  ■ 


The  homogeneous  forms  of  (4.15a,b)  have  the  eigeiisolutions 


Uo  =  Cl  ^  exp 


9 


Oo  =  C2^  exp 


(4.16) 


where  Ci,C2  are  arbitrary  constants.  It  follows  that,  since  the  inhomogeneous  solutions 
of  (4.15)  must  satisfy  Uo  =  0,  =  0  at  ^  =  0,  the  algebraically  decaying  solutions  of  the 

homogeneous  forms  of  (4.15a,b)  are  recpiired.  However  if  these  solutions  are  retained  we 
cannot  match  with  the  corcflow  so  t>o  nmst  be  zero.  Hence  the  highest  order  term  in  the 
expansion  of  v  is  vi  and  (2.7)  then  gives 


(4.17) 


-A' 


which  has  the  solution  vy  =  —^^Uo  satisfying  the  boundary  conditons.  The  functions 
Uo,  fi,  &o  are  shown  in  Figure  (4.1).  At  leading  order  in  the  expansions  of  it,  v  and  0  it 
is  clear  that  the  disturbed  flow  is  confined  to  the  wall  layer  region,  however  at  next  order 
the  flow  is  no  longer  contained  in  the  wall  layer.  The  function  wi  satisfies  the  equation 


“■  +  y"‘  +  3  ‘  " 


(4.1Sa) 


and  hence 


-3Cin'^  [  ui(0 


where  rii  is  the  exponentially  decaying  solution  of 


(O)J’ 


(4.186) 


£ 

-  ,  s  -  /  ,  s  >4  n 

U]  +  — »i  +  -«i  —  0, 


IS 


and  6i  is  similarly  given  by 


-3020^ 


^i(0] 

^i(0)i 


Here  Oi  is  the  exponentially  clecH3’ing  solution  of 


1  -  f  - 

(7  0  0 


(4.19) 


At  higher  order  ?’2  is  found  to  satisfy 


4  +  +  Y  )  '’2  =  +  V'’* 


2^ 


2A' 


3A 


Equation  (4.20)  is  solved  subject  to  the  boundary  conditions  e2  =  i’2  =  0 ,  C  =  0 
^  =  constant  +  exponentially'  small  terms  as  — >  (X).  Integration  of  (4.20)  numcTically 
shows  that  this  constant  is  nonzero,  hcnct'  the  perturbed  velocity  and  temperature  fields 
for  the  y  =  0(  1 )  region  arc  (jf  the  form 


u  ^  --  j„{  y,  a)  ,  r~7??(j/,a),  (4.21) 

a  ^ 

where  m  satisfies  (4.1 1 ).  Hence  it  follows  that  ( </,  v,0)  >  0  as  y  ►  oo.  The  function  Tn{y,  a) 

is  shown  in  Figure  (4.2)  for  diffen-nt  vortex  wavenumbers.  By  combining  (4.13,4.20)  a 
composite  disturbance  field  is  obtained  for  .some  small  value  of  x  which  can  be  used  as 
the  initial  condition  for  the  solution  of  the  full  linear  disturbance  equations  (2.66),  (2.6e), 
(2.7),  subject  to  the  boundary  conditions  (2.7,/ ). 

The  calculations  we  nf)w  reprwt  on  were  carried  out  u^ing  the  numerical  scheme  de¬ 
scribed  earlier  with  the  starting  point  of  the  calculation  at  x  =  0.505  and  with  the  stepsize 
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in  the  x  direction  e  equal  to  0.004.  In  order  to  determine  the  effects  of  the  wall  forc¬ 
ing  functions,  imposed  at  the  position  x  =  on  the  flow  downstream  of  that  point  the 
dimensionless  energy  E{x)  of  the  flow  was  monitored  where 

rOO 

E=  {u^{x,ij)  +  u^{x,ij)  +  w'^(x,y)}  dy 
Jo 

Q  E  j  Q  V 

and  the  local  growth  rate  o{x)  —  = —  .  The  position  of  neutral  stability  is  defined  as  the 

E 

position  at  which  a  =  0.  Of  course  other  instability  criteria  can  be  defined  but  some  limited 
experimentation  showed  that  the  neutral  curve  is  not  greatly  dependent  on  the  choice  of 
flow  property  used  to  monitor  the  growth  of  the  disturbance.  Moreover,  we  believe  that 
the  flow  property  we  have  used  is  a  sensible  one  because  it  accounts  for  the  changes  in  all 
of  the  velocity  components  in  some  sense  averivged  across  the  flowfield.  For  a  given  wall 
forcing  function,  G  and  different  values  of  a  we  marched  downstream  and  calculated  the 
position  at  which  the  vortex  structure  began  to  grow.  The  local  wavenumber  and  the 
local  Grashof  number  Gj  were  calculated  and  a  neutral  curve  in  {aj:,Gi)  space  formed. 
This  process  was  repeated  for  different  values  of  G.  Two  sets  of  initial  conditions  were 
considered.  Firstly  the  problem  was  solved  for  an  initial  disturbance  of  the  form 

u  =  ^  ~  ’  0  —  Q  imposed  at  x  =  x  -f  .005.  (4.22) 

This  type  of  disturbance  correspondes  to  the  situation  when  the  vortices  are  stimulated  by 
wall  roughness,  later  we  shall  look  at  the  case  when  the  vortices  are  induced  by  non-uniform 
wall  heating. 

The  development  of  the  velocity  and  temperature  fields  downstream  is  shown  in  Fig¬ 
ures  (4.3a-c)  for  G  =  5,  a  —  0.4.  The  corresponding  neutral  curves  for  the  problem  are 
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shown  in  Figure  (4.4).  We  see  that  as  G  varies  the  neutral  curves  move  around  in  the 
wavenumber- Grashof  number  plane.  Our  calculations  show  that  there  is  a  neutral  curve 
corresponding  to  G  near  2  which  has  the  lowest  minimum.  Furthermore  the  minimum  value 
of  the  Grashof  numljer  on  that  curve  is  about  4.  When  the  Grashof  number  is  increased 
or  decreased  from  this  value  the  curves  move  upwards  so  the  the  flow  is  more  stable.  It 
should  be  noted  that  the  effect  of  changing  G  in  the  calculation  of  the  neutral  curves  shown 
in  Figure  (4.4)  is  exactly  equivalent  to  calculations  carried  out  with  a  fixed  G  but  with 
the  position  where  the  forcing  Ix'gins  now  being  varied.  Thus  in  Figure  (4.4)  increasing  G 
corresponds  to  moving  the  forcing  further  downstream  from  the  leading  edge.  In  that  case 
it  is  not  surprising  that  the  curves  in  the  Oi  —  Gi  plane  move  upwards  since  the  forcing 
cannot  initially  generate  an  unstable  vortex  since  the  form  of  the  initial  disturbance  is  not 
typical  of  a  growing  vortex  flow.  Likewise  when  G  is  decreased  in  Figure  (4.4)  the  forcing 
is  being  moved  progressively  towards  the  leading  edge,  the  fact  that  the  curves  move  up¬ 
wards  in  this  case  again  implies  that  rt)ughuess  near  the  c?dgc  has  a  relatively  weak  effect 
on  the  flow.  In  an  experiment  one  would  exp<'ct  that  localized  roughness  elements  would 
be  distributed  at  several  sites  along  the  wall  so  that  the  most  dangerous  mode  would  be 
the  one  excited.  We  postpone  a  discussion  of  the  available  experimental  results  until  the 
final  section  of  this  piiper. 

In  fact  the  discussion  given  above  applies  equally  well  to  the  case  when  the  forcing 
corresponds  to  a  localized  temperature  \'ariation  at  the  plate.  In  particidar  for  the  problem 
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with 


imposed  at  a;  =  x  +  .005. 

(4.23) 

We  now  find  that  the  most  dangerous  disturbance  occurs  when  G  is  close  to  8  whilst 
the  local  Grashof  number  corresponding  to  this  most  dangerous  mode  is  about  2.  The 
downstream  development  of  the  velocity  and  temperature  fields  is  shown  in  Figures  (4.5a- 
c).  The  corresponding  neutral  curves  are  shown  in  Figure  (4.6).  Again  for  large  values  of  G 
we  can  see  that  the  forcing  applied  docs  not  initially  generate  unstable  vortex  structures.  It 
is  also  clear  that  the  forcing  becomes  less  <langerous  when  it  is  moved  close  to  the  leading 
edge.  In  each  of  the  above  calculations  it  was  found  that  the  dimensionless  energy  E 
decreases  by  several  orders  of  magnitude  before  the  growth  rate  a  becomes  positive.  This 
means  that  the  forcing  applied  generates  vortices  which  decay  significantly  before  they 
begin  to  grow.  Hence  a  localized  wall  forcing  function  is  not  a  particularly  efficient  means 
for  the  production  of  longitudinal  vortices.  However  note  that, despite  being  an  inefficient 
generator  of  longitudinal  vortex  rolls,  in  the  absence  of  other  forcing  modes  an  isolated 
wall  heating  forcing  will  result  in  vortex  growth  downstream.  We  conclude  this  section  by 
noting  that  if  instability  is  caused  by  isolated  wall  roughness  or  nonuniform  wall  heating 
then  we  expect  instability  will  occur  when  the  local  Grashof  number  exceeds  about  2. 

§5  Free-stream  Disturbances 

We  shall  now  consider  the  generation  of  vortex  structures  due  to  a  freestream  lon¬ 
gitudinal  vortex  field  imi)inging  on  the  h'ading  edge  rather  than  imposing  some  initial 


=  0  , 


V  — 


0  =  ^exp 
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disturbance  on  the  flow.  We  take  the  streainwise  velocity  component  to  be  of  the  form 


(/  =  1  +  A(:“‘~Uc(u)  ,.r  =  0, 


(5.1) 


and  hence  we  have  assumed  a  dependcnc<*  of  the  impinging  vortex  structure  on  the  bound¬ 
ary  layer  lengthscale.  However  it  will  later  be  shown  that  the  case  of  Uc  independent 
of  y  gives  rise  to  the  most  dangerous  vortfcx  mode.  We  consider  Uc  to  be  of  the  form 
Uc{y)  =  cos(by  +  (i>)  where  b  and  4>  ar<^  constants  so  that  the  disturbance  is  periodic  in  both 
the  y  and  2  directions.  We  need  to  consider  two  regions,  the  boundary  layer  y  ~  and 
also  an  outer  region  y  ~  0(1),  this  is  because  the  wavelengths  in  the  spanwise  and  normal 
directions  are  large  compared  with  the  boundary  layer  scale  ,t^. 

Consider  the  boundary  layer  region  y  ~  .r^  with  .c  <<  1  and  allow  y/\/x  — >  oo,  the 
disturbance  ecpiations  (2.G,  2.7)  may  be  written  in  the  form 


u  =  0. 
^  =  0, 


where 

d  =  lim  ()//'  -  /), 

— OO 

is  the  Blasivis  constant. 


(5.2a) 

(5.26) 

(5.2c) 

(5.3) 


We  require  a  solution  of  (5.2a-c)  which  will  maintain  the  periodicitv  —  in  the  y 

0 

direction,  hence  we  take 


(5.4a) 


u  =  u„  =  c  "  cos(bij  +  o  —  b;i\/'2.r). 


$  =  6p  =  O^c  cos(/a/  +  0  —  b,i'i\^^‘). 


a-GOia 


c-i{l  -  a) 


cos{hy  +  0  —  /;.i\/2r)  for  a  ^  1, 


[5.4c) 


I  a-  sj2x 

e~“  ■'  cos{htj  +  (f)  —  bd\/^-)  for  <7  =  1, 


(5.4f/) 


where  B\,k\  and  k-i  are  arbitrary  con.stants. 


There  is  an  eigensolution  of  (5.2c)  ii  =  Q(x}e  for  arbitrary  Q{x)  and  this  solution 
is  needed  in  order  to  match  with  the  lioundary  layer  solution.  Hence  the  appropriate 


solution  of  (5.2c)  is 


V  =  V,,  +  Q(x)e' 


and  it  is  to  be  noted  that  the  Iroundary  layer  .structure  causes  the  periodic  form  of  r  to 
occur  only  when  ay  >>  1.  For  the  region  where  .r  =  0(1)  and  a?/  >>  1,  at  the  edge  of  the 
boundary  layer,  we  must  solve  the  disturl)auce  etjuations  subject  to 


[a.  r. H)  =  [Up.  Vp.Bp] 


\\v  now  considf'r  the  ;/  =  0(1)  region  with  x  «  1  and  w('  (h'tc'rmine  Q[x)  by  matching 
with  the  boundary  layer  solution.  Tlu'  functions  u.a  and  B  arc'  obtained  from  expansions 
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in  powers  of  xi  in  the  disturbance  cfpiations  by  perturbing  the  basic  flow  in  the  form 


u  =  cos(^{/'  +  y/2f"]  +  •  •  • , 

(5.7a) 

B  =  cos  0{/j'  +  ij/'I}'."}  +  ■  •  ■ , 

(5.76) 

(5.7c) 

where  /?(>/)  satisfies  Ji'"  ~  —afh". 

(5.7d) 

and  as  ?7  — +  oo 


cos  0/9 


Matching  with  the  boundary  layer  solution  for  rj  >>  1  yields 

(^9  -  6/9^)  [f  (6^  -  a^)  .  '  a-Gdia 

Q  +  cos - 5 -  - 7k=~^  ■'  '  -  r  +  ~Y? - TT^ 


13  cos  4> 

~  2v/2x 

A  composite  solution  is  formed,  from  the  y  ~  x^  and  y  ~  0(1)  solutions,  for  small  x  to 
give  the  asymptotic  forms  for  u,e  and  B,  and  hence  initial  data,  for  the  solution  of  the 
disturbance  ecpiations  using  the  numerical  scheme  starting  at  a  small  value  of  x. 

We  note  that  when  (;i»  ==  0.  the  function  a,  takes  the  form  —  cos{by)  which  means 
that  the  incoming  vortex  field  does  not  satisfy  the  iio-slip  condition  at  tlic  wall.  Hence  we 
choose  4>  to  be  zero  since  this  relates  to  the  most  physically  relevant  case  corresponding  to 
Uc(0)  ^  0.  We  again  monitored  the  dimensionless  niergy  £'(x)  of  the  flow  where 


E  =  f  {a^(x,  y)  +  e^(x,  y)  +  a'^(x,  y)}fly. 
Jo 


QE I  ^  V 

and  the  local  growth  rate  <7(x)  =  - —  .  The  j)osition  of  neutral  stability  is  ch'flned  as 

E 

the  position  at  which  a  =  0. 


For  the  receptivity  problem  fornnilated  above  the  disturbance  equations  were  marched 
downstream  using  the  numerical  sclunne  described  in  Section  3.  The  steplength  in  the 
streamwise  direction  was  taken  to  be  0.00001,  this  very  small  x  steplength  was  necessary 
due  to  the  singular  behaviour  of  u  for  small  x. 

The  profiles  for  the  disturbance  velocity  and  temperature  components  as  the  vortex 
develops  downstream  are  shown  in  Figures  (5.1a-c).  We  can  see  that  the  edge  velocity  for 
the  streamwise  velocity  component  decreases  monotonically  with  x  as  does  the  temperature 
component.  The  normal  velocity  component  v  at  the  edge  of  the  boundary  layer  is  seen 
to  increase  as  we  move  downstream, however  for  larger  values  of  x  the  edge  velocity  begins 
to  decrease  with  x  due  to  the  exponential  factor  in  (5.4c).  In  Figures  (5.2)-(5.3)  the  local 
Grashof  number  and  the  local  wavenumber  have  been  calculated  at  the  points  of  neutral 
stability  where  the  local  growth  rate  vanishes;tliis  enables  us  to  generate  neutral  curves  in 
the  local  wavenumber  -  Grashof  number  space.  For  these  two  cases  the  parameters  chosen 
were  G=70,^i  =  l.  We  conclude  from  these  calculations  that  instability  first  occurs  when 
b=0,  in  fact  further  calculations  for  different  values  of  h  produced  neutral  curves  located 
above  that  for  b  =  0.  Wc  deduce  that  Ur  ~  cosaz  is  the  most  dangerous  form  for  the 
incoming  vortex  field.  We  conclude  from  Figure  (5.2)  that  freestream  disturbances  are 
able  to  cause  the  onset  of  instability  when  the  local  Grashof  number  is  bigger  than  about 
.02;  this  is  significantly  lower  than  the  critical  Grashof  number  associated  with  isolated 
wall  forcing. 

Distributed  Roughness 


We  shall  now  consi^'  j.  ^he  effects  of  wall  heating  on  the  forced  convection  flow  when 
the  wall  forcing  is  described  by  a  non-localized  forcing  function  and  occurs  on  an  0(1) 
strcamwise  lengthscale.  Again  we  concentrate  on  the  0(1)  wavenumber  regime.  The  lin¬ 
earized  disturbance  equations  (2.G6),  (2.6e),  (2.7)  were  solved  together  with  the  boundary 
conditions  (2.G/)  where  Fi{x)  is  given  by 

Fi(.t)  =  40 -  i^,.-20(x-i/2)=  ^ 

and  jF2  =  0  so  that  the  vortex  is  induced  by  wall  loughiress  rather  than  non-uniform  wall 
heating.  Given  a  function  q{a)  the  disturbance  ccjuations  can  be  marched  downstream 
from  X  =  0  for  various  values  of  a  using  the  numerical  scheme  described  in  section  3  with 
no  initial  disturbance.  We  then  invert  this  transformed  flowfield  in  {a,x,y)  space  to  give 
the  flowfield  in  {z,x,y)  s^iace  induced  by  a  forcing  function  of  lieight  proportional  to  q{z) 
the  inverse  Fourier  transform  of  q{a).  A  symmetric  obstacle  q(z)  was  considered  with 

<;(--)  =  ^cx,>(^-^j,  (6.2) 

and  q{a)  the  Fourier  transform  of  (G.2)  is  given  by 

=  ^exp(-4o^).  (6.3) 

We  then  combine  (G.l)  and  (C.3)  to  give  the  boundary  conditions  (2.C/)  in  (a,x,  y)  space. 
The  disturl)ance  eciuations  were  marched  downstream  and  t?*(r/,.r,y),  the  maximum  value 
of  9,  was  calculated.  This  procedure  was  repeatc'd  for  various  values  of  a  and  the  transform 
in  r  was  then  inverted  numerically  to  give  9{z,x,  y).  The  parameters  chosen  were  x  varying 
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between  0  and  10  whilst  ^  varies  between  0  and  10.  The  calculation  was  repeated  for 
different  values  of  the  Grashof  number  G. 

The  velocity  and  temperature  profiles  for  the  syiiimetric  obstacle  are  shown  in  Figures 
(6.1a-c)  for  G  =  S,  a  =  0.45.  It  is  to  be  noted  that  the  nia.x'imum  value  of  d  occurs  initially 
at  the  wall  and  as  x  increases  its  position  moves  away  from  the  wall. 

For  a  distributed  wall  forcing  function  there  is  a  strong  coupling  between  the  induced 
vortex  field  and  the  wall  forcing  with  the  ratio  between  the  two  being  a  function  of  the 
wavenumber  and  the  Gbrtler  number.  In  principle  we  could  maximise  the  coupling  between 
the  vortex  field  and  the  forcing  by  varying  a  and  G  though  this  would  require  a  large 
amount  of  computer  time. 

The  contours  of  9*{z,x,y)  arc  shown  in  Figure  (6-2).  They  demonstrate  that 
immediately  after  the  obstacle  the  perturbed  temperature  field  decays  and  is  formed  into 
a  wake  solution  behind  the  obstacle.  However,  further  downstream  the  effects  of  thermal 
instability  due  to  the  heated  wall  reainplifics  this  disturbance  into  longitudinal  vortex  rolls. 
The  distance  between  the  reamplification  of  the  clisturbance  and  the  ol^stacle  decreases 
as  we  increase  the  Grashof  number.  The  same  effect  could  have  been  demonstrated  by 
following  the  same  method  as  in  section  5  by  fixing  the  Grashof  number  and  varying  the 
position  at  which  the  forcing  was  first  applied. 

The  ?»me  calculations  were  repeated  for  an  asymmetric  oljstacle  of  the  form 

=  .  (0^4) 

The  H'sults  are  shown  in  Figure  (6.3)  .  A  similar  flow  structure  to  that  of  the  symmetric 
obstacle  was  obs(;rved  with  the  perturix'd  flowfi(dd  formed  into  a  wake  l)efore  subsetpiently 


2S 


being  rcamplified  further  downstream  into  longtitudinal  vortex  roll  structui'es.  Again 
the  distance  of  the  resulting  vortex  structure  from  the  leading  edge  is  dependent  upon 
the  Grasliof  number.lt  is  worth  noting  that  similar  flow  structures  have  been  reproduced 
experimentally  by  Mangalani  et  al.(19S7)  for  the  analogous  Gdrtler  type  vortex  problem. 
§7  The  high  Grasliof  number  limit  and  the  fastest  growing  mode. 

At  high  values  of  tlie  Grashof  number  we  expect  that  viscous  effects  will  be  negligible 
except  at  low  or  high  vortex  wavenumbers.  An  examination  of  (2.6)  for  G  >>  1  with  a 
held  fixed  suggests  that  inviscid  disturbances  have  ~  G'^;  we  are  therefore  led  to  the 
expansions 


[u,G  2u,G  '^iv.G  *i>,  0]  =  [(wo^  t’o,  u’o.po,  ^0 )  +  •  • C^i3{z)di^ 


(7.1) 


where  etc.  are  functions  of  x  and  y  only  whilst  3  expands  as 


3  —  3o  3-  3iG  2  4- 


(7.2) 


If  the  above  expansions  are  substituted  into  (2.G)  and  the  dominant  terms  are  retained  in 
the  limit  G  ^  oo  we  obtain 

do'f'o  +  i'ojj  +  idu'o  =  0, 

3(}II(Id  +  I'olly  =  0. 

=  —poy  +  ^0.  (7.3n  —  e) 

3()d9[)  4-  <'\)Ty  =  0, 
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and  we  can  eliminate  «0i  J'lul  fiom  this  system  to  give 


u[v{)yy  -  o‘(’oj  -  Uyyi'o  =  (7.4) 

PqU 

This  equation  must  be  solved  subject  to  i>o  =  0  at  y  =  0,  oo  and  this  specifies  an  eigenvalue 
problem  do  =  du(«). 

Thus  in  the  inviscid  limit  the  growth  of  the  disturbance  is  governed  by  a  quasi-parallel 

stability  problem  since  in  that  case  the  disturbance  varies  on  a  relativ'ely  short,  G~K 

lengthscale  in  the  x  direction.  We  shall  restrict,  our  attention  here  to  the  determination  of 

do;  higher  order  terms  in  the  <^xpa.nsion  of  the  growth  rate  can  be  obtained  in  a  routine 

manner.  We  note  that  we  have  assumed  that  there  are  unstable  solutions  of  (7.4),  this  is 

assured  if  there  are  regions  where  Tj,  <  0  in  0  <  y  <  oo.  The  numerical  solution  of  the 

eigenvalue  problem  specified  by  (7.4)  together  with  the  conditions  uo  =  0,  y  =  0,  oo  is  made 

nontrivial  because  of  the  singularity  in  the  equation  at  y  =  0.  An  examination  of  (7.4) 

for  y  <<  1  shows  that  for  y  <<  l,eo  ~  y^.  In  fact  (7.4)  is  more  easily  solved  by  making 

the  transformation  ;/  =  — — — —  log  y  and  the  results  presented  below  were  obtained  using 

Uy(.r,0) 

that  transformation. 

In  Figure  (7.1)  we  show  the  most  unstable  eigenvalue  for  <7  =  1,  it  can  be  seen  that 
the  growth  rate  increases  monotonically  with  a.  The  eigenvalue  shown  is  in  fact  just  one 
of  an  infinite  secpience  of  unstable  modes.  For  small  value  a  we  see  that  the  growth  rate 
goes  to  zero  like  some  powers  of  a;  actually  our  calculations  suggest  that  do  ~  for 
a  <<  1  and  this  asymptotic  limit  will  be  considered  later  in  this  section.  In  Figure  (7.2) 
we  have  shown  eo(y)  for  three  differc’iit  values  of  a;  we  note  that  the  disturbance  becomes 


less  ccjiicentrated  as  a  decreases.  Now  we  investigate  further  tlu*  inviscid  problem  at  large 
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wavenvvmbci's.  This  will  euabh'  \is  to  identify  the  scale  on  which  viscons  effects  come  into 
play  and  thci'(Tore  we  will  be  in  a  position  to  identify  the  fastest  growing  mode. 

For  large  values  of  a  the  effect  of  normal  diffusion  of  heat  and  vorticity  can  only  be 
comparable  with  that  in  the  spanwis(>  direction  if  the  vertical  depth  of  vortex  activity  is 
compressed  to  0(u“'  ).  We  shall  therefore  seek  a  solution  of  (7.4)  within  a  hu'er  of  depth 
at  the  boundary:  we  therefore  dc'fine 


(7.5) 

and  note  that  n,T  expand  as 

u  =  nC,a~^  +  •  •  • ,  T  =  1  —  tc(^n“*  +  •  •  •  .  (7.6a,  6) 

In  order  that  the  dominant  terms  on  the  left  hand  side  of  (7.4)  are  comparable  with  the 
term  on  the  right  hand  side  we  must  write 

^0  =  U/doo  +  •  •  •  •  (7-7) 

The  zeroth  order  approximation  to  the  eigenvalue  problem  for  /3o  then  reduces  to 

yi  2 

-  '’u}  =  -coC~".  co  =0,  C  -  0,oo.  (7.S) 

The  solution  of  the  above  equation  can  be  e.xjnessed  in  terms  of  modified  Bessel  functions 
of  imaginary  orcU’r.  How('ver  we  were*  unabh'  to  find  the  recpiired  asymptotic  properties  of 
such  functions  so  we  must  solve  it  numerically.  Since  out  i)rimary  aim  at  this  stage  is  to 
see  how  viscous  effects  come  into  play  it  is  not  ueci'ssary  for  us  to  solve  (7.S);  in  fact  we 
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can,  if  necessary,  infer  the  vahu'  of  /ioo  from  the  limiting  small  wavenumber  approximation 
to  the  viscous  calculations  to  be  discussed  next. 

The  large  a  inviscid  analysis  given  above  must,  of  course,  breakdown  when  iidx  r\j  dy^ 
in  which  case  viscous  effc’cts  cannot  be  ignored.  This  balance  is  achieved  where  a  ~  G  '^  so 


we  now  write 


(I  =  G  ‘  a , 


and  modify  (7.1)  to  give 


[a,G  ^v,G  ^w,G  •' p,  d]  —  [(  i/q,  t-oed’o,  75o,  ^0 )  +  •  •  ^ 


where  etc.  are  functions  only  of  .r  and  (,'  defined  by  (7.5).  The  eigenvalue  /3  now 

expands  as 


a  =  1^(1  +  iiiG  •*  + 


(7.10) 


and  the  zeroth  order  growth  rate  is  found  to  be  determined  by  the  eigenvalue  problem 

^odo  T  dcoc  +  iaii'o  —  0. 


r 

^'0  =  Tyt’o, 

_  1  ,  ^0 
^e()  —  —  — 

(I  a- 


^>i'o  —  — Po¬ 
et 


ra  - 

—  r.,'  <  ()• 


>io  =  C(|  =  (Co  =  Ou  =  0,  C  =  (hoc- 


Her(’  the  o])erator  C  is  defiiu'd  by 


(7.11e-/) 


L  =  0^  -  1  -  ~,io- 

(I  * 
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A  more  convenient  form  of  the  disturbance  equations  can  l^e  obtained  by  eliminating  po 
and  Wo  to  give 


e  ^  \ A  (yw  * 

Cuq  —  -rr-,  —  l)i’o  =  TV,  Cuo  —  t'o 


uo  =  Vo 


Vo  =  ^0  —  0,  C  =  Ih  oc 


(7.12a  -  d) 


A  point  which  should  be  noticed  here  is  that  the  i  momentum  equation  decouples  from  the 
other  equations  so  that  the  eigcnrelation  is  determined  by  the  sixth  order  system  associated 
with  Vo,  do-  It  is  also  possible  to  scale  /i,a  and  iv  out  of  the  above  eigenvalue  problem  by 
redefining  d  and  /3.  In  Figure  (7.3)  we  show  the  most  unstable  eigenvalues  of  (7.12);  the 
results  shown  were  obtained  using  a  fourth  order  finite  difference  scheme  to  discretize  the 
differential  equations  for  vq  ‘ind  do. 

For  small  values  of  d  we  see  that  /So  ~  d  so  that  we  obtain  the  required  match  with  the 
large  wavenumber  limiting  form  of  the  inviscid  mode.  The  growth  rate  attains  a  maximum 
at  a  finite  Vcilue  of  d  and  then  passes  through  zero  at  a  sufficiently  large  value  of  d.  In  fact 
this  zero  of  the  growth  rate  corresponds  to  the  right  hand  branch  of  the  neutral  curve  in 
the  Grashof  number  -  wavenumber  plane.  Actually  it  is  only  possible  to  find  solutions  of 
(7.12)  with  fSo  7^  0;  the  case  (So  =  0  corresponds  to  the  ca.se  when 

-  _  d  ~  aw  ^ 

—  Vo  —  —7,  —VO  —  TrTTCo 
a-  cw 

so  that  d"*  =  ato  which  means  that  the  zeroth  order  approximation  to  the  right  hand  branch 
of  the  neutral  curve  is  given  l;y 

G  =  —  +  .  (7.13) 

(7UJ 
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An  investigation  of  (7.12)  with  Bq  <<  1  shows  that  the  wall  layer  becomes  thicker  in 
this  limit  and  higher  order  terms  in  (7.13)  can  be  found  when  it  has  become  0(G5);  the 
structure  in  that  case  is  similar  to  that  found  by  Hall  (19S2a)  for  the  Taylor  problem. 

In  Figure  (7.4)  the  eigenfunctions  me  shown  for  three  different  values  of  the 

wavenumber.  We  note  that  at  the  largest  value  of  the  wavenumber  the  eigenfunctions  have 
spread  further  away  from  the  wall  whilst  at  the  smallest  wavenumber  the  temperature 
disturbance  develops  a  wall  layer  structure  consistent  with  the  inviscid  limit  discussed 
above.  Thus  we  have  shown  above  that  at  high  Grashof  numbers  the  fastest  growing 
disturbance  is  localized  at  the  wall  and  is  dominated  by  viscous  effects.  The  unstable  band 
of  wavenumbers  cuts  out  at  the  right  hand  branch  of  the  neutral  curve  where  a  ~  (aa-'G)^. 
For  a  ~  0(1)  the  disturbances  are  essentially  inviscid  with  a  growth  rate  which  leads  to  zero 
where  a  — >  0.  At  some  stage  viscous  effects  will  reappear  for  sufficiently  small  a  and  then 
we  expect  that  the  left  hand  branch  of  the  neutral  curve  will  be  encountered.  Though  the 
growth  rates  in  this  regime  are  relatively  small  it  is  important  for  disturbances  localised 
very  close  to  the  leading  edge  of  the  wall.  In  addition  there  is  an  unexpected  connection 
between  vortex  disturbances  and  Tollmien-Schlichting  waves  here  so  it  is  now  considered 
in  some  detail.  As  a  first  step  we  consider  the  limiting  form  of  (7.3)  where  a  — +  0.  It  is 
clear  from  (7.4)  that  we  must  consider  .separately  the  regions  y  =  0(1)  and  y  =  0(a~^). 
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For  2/  =  0(1)  we  write 


XL  —  XIq  it  Hi  + 


V  =  i?o  +  rt  2  i?j  -|- 


i  ^ 

w  =  a^WQ  +  a^wi  +  '  • '  , 


(7.14a  -  e) 


and  /?o  then  expands  as 


^  —  ^0  +  0.^1  +  •  •  • , 


p  =  2^0  +  «Pi  +  •  •  • , 


+  /d'a  +  •  -  •  . 


(7.15) 


If  we  now  substitute  the  above  expansions  into  (7.3)  and  solve  the  leading  order  approxi¬ 


mation  to  this  system  'vr  Ind  that 


Uq  —  ^ y 1 


ro  = 


^0  =  -r„ 


(7.16a  —  d) 


Po  =  -T  +1. 

Meanwhile  in  the  upper,  0(a~* ),  layer  we  can  easily  show  that 

p  =  ••• 

V  =  - . 

Thus  we  can  only  match  the  expansions  for  a  in  the  two  layers  if 


(7.17a,  b) 


(/?")'  =  1 


so  that  for  small  a  the  eigenrclation  associated  with  (7.3)  takes  the  form 


/^o  =  (1^  + 


(7.18) 
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The  small  a  inviscicl  solution  discussed  above  fails  when  udx  ~  in  the  viscous  layer  of 
depth  a~^  adjacent  to  y  =  0.  This  occurs  when  a  ~  G'“?  and  we  then  have  a  triple-deck 
like  structure  with  three  layers  of  depth  OICt"" ^  ),  0(  1 ),  0(G' ^ )  to  consider.  (Note  that  the 
structure  described  below  applies  also  to  the  Gbrtler  problem  and  was  alluded  to  by  Denier, 
Hall,  and  Seddougui  (1991).)  We  therefore  write  a  =  G~"^a  . 

In  order  to  allow  for  the  possibility  of  unstable  Tollmien-Schlichting  waves  we  must 
modify  (2.C)  to  allow  for  the  possibility  of  time-dependent  modes.  This  is  simply  done  by 
inserting  the  term  u_i  into  the  momentum  efiuations  and  the  term  Bt  into  the  temperature 
equation.  In  the  lower  deck  where  y  =  0(G'“^ )  we  define  ^  and  write 

{u,G~^  iuG~hv.G~^  p,B}  ==  {(t/o,t’o,n’o,Po,^o)  + 

where  jj  =  ^qG^  -f  •  •  • ,  Q  =  flG?  -!-•••,  whilst  uo,  t’o,ctc.  are  functions  oi  ^,x.  We  assume 
that  the  frequency  Q  of  the  disturbance  is  constant.  The  equations  to  determine  the  zeroth 
order  approximation  to  the  disturbance  in  the  lower  deck  are  found  to  be 

+  t’Qc  +  ifiu'o  =  0, 

—  iSliiQ  -j-  00(1^110  T  t’o/i  = 

0  —  po‘,  (7.19a  — c) 

-Hhc()  -t-  /Aj//G'’o  =  ->Po(>  + 

—  iQBo  -f-  t%PsBo  +  *-’0^’  = 

(7 

which  must  be  solved  subject  to 

(/„  =  )'(,  ^  (Co  =  B^f  =0,  (f  =  0,  (7.20) 
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and  appropriate  matching  conditions  at  ^  =  oo.  Tims  in  the  lower  deck  we  have  that 

P  =  -Po, 

— (i"Po 

+  i(in'i)  =  — ^ ; -  /  Ai{^)d(p, 

A3T,(0o)  ^00 

rvhcre  Pq  is  a  coustaut,  (f)  —  =  <P{^  =  0),  A  =  f^op  and  .4,  is  the  Airy 

function.  In  the  main  and  upper  decks  the  disturbance  takes  on  essentially  the  same  form 
as  that  discussed  above  for  the  small  a  limit  of  the  inviscid  problem.  In  particular  we  find 


that  in  the  main  deck 


a 

u  =  —Cuy, 


where  C  is  a  constant. 


Thus  matching  l^etween  the  main  and  lower  decks  is  achieved  if 


c  [4  - 1] = po,  = -  r~ r 


which  leads  to  the  eigenrclatioii 


[-1  +  /  A.df]  =  A5A'( 


(7.21) 


We  note  that  if  we  take  the  further  limit  a  oc  we  recover  the  limiting  inviscid  solution 
(7.18)  whilst  in  the  limit  a  — >  0  we  obtain 


S'^a  f  .4,(/// = 
J  (f>(\ 


(7.22) 


and  this  eigenrelation  can  be  found  from  the  limiting  large  wavenumber  analysis  of  Hall 


and  Smith  (1984).  In  fact,  rather  than  .solve  (7.21)  for  .do  as  a  function  of  a,  it  is  more 
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instructive  to  rewrite  this  equation  using  the  inverse  of  the  spanvvise  wavenumber  rather 
than  G  as  the  appropriate  large  parameter.  In  order  to  do  this  we  write 


in  which  case  (7.21)  becomes 

+  =  (7.23) 

where 

/OC 

^  A,(,?)di/,.4:,  =  (7.24) 

If  we  set  G'o  =  0  in  (7.23)  the  resulting  equation  determines  the  scaled  growth  rate  of  a  3D 
Tollmien-Schlichting  wave  of  freciucncy  0,* .  By  varying  G'o  in  the  range  —  oo  <  G'o  <  oo 
we  can  then  infer  the  effect  of  heating  or  cooling  on  very  oblique  Tollmien-Schlichting 
waves.  Alternatively  by  setting  Q*  =  0  we  can  obtain  the  required  match  with  the  inviscicl 
low  wavenumber  modes  discussed  ('arlier.  In  fact  we  can  see  directly  from  (7.24)  that,  if 
Go  >>  1  with  Q*  held  fixed,  then  the  zeroth  order  approximation  to  the  eigcnrelation  is 
simply 

i3q  =  ( Go )  ^  +  •  •  • ,  Go  >  0.  ( 7.25) 

This  corresponds  to  the  limiting  inviscid  form  (7. IS)  and  we  also  deduce  that  instability 

occurs  only  for  positive:  Grashof  numbers.  Now  we  shall  present  results  for  the  solution  of 

(7.23)  for  a  range  of  values  of  G.  It  is  well-known  that  neutral  solutions  of  (7.23)  occur 
A'  i 

when  — —  =  Ni^'  with  N  ~  1.001.  Tims  the  neutral  values  of  3*  are  given  by  3*  =  7  3\, 


The  solutions  of  this  equation  are  shown  in  Figure  (7.5).  The  frequency  Qn  associated 


with  each  of  the  neutral  points  is  given  by 

~  -/h.29S/}|. 


Thus  neutral  solutions  exist  only  for  nonzero  frequencies  although  we  note  that  for  Gq  < 
0,  |Gu|  >>  1  equation  (7.23)  yields 


d*  =  ?lGol^  + 


.  i 

t  a 


-1 


The  leading  order  term  here  reproduces  the  leading  order  neutral  inviscid  result  appropriate 
to  negative  Grashof  numbers  whilst  the  second  viscous  term  always  has  negative  real  part 
so  that  the  flow  is  stable  in  that  case.  The  latter  results  holds  for  all  frequencies. 

In  Figure  (7.0)  we  show  the  dependence  of  /3*  on  the  frequency  12*  for  a  range  of 
values  of  the  scaled  Grashof  number  G'o-  We  see  that  the  mode  with  Go  =  0  is  unstable 
for  12*  >  2.298,^*  >  1.  Between  12*  =  0,2.298  the  mode  is  stable  but  we  note  that  the 
growth  rate  approaches  zero  when  12*  —*  0.  If  Go  is  now  taken  to  be  slightly  positive 
then  this  limiting  neutral  point  moves  to  a  small  positive  value  of  12*  and  there  is  then  a 
small  but  finite  band  of  unstable  distuibances  of  small  frequency.  When  Go  is  increased 
beyond  about  .00  this  unstable  band  connects  with  the  other  unstable  band  corresponding 
to  the  modes  with  12*  >  2.298  for  G  =  0.  At  this  stage  there  are  no  neutral  solutions 
and  as  Go  is  increased  the  growth  rates  at  small  frequencies  increase  and  approach  the 
limiting  case  0*  =  IGol^  for  Go  >>  1.  However,  at  any  fixed  value  of  Go,  we  see  that  the 
growth  rate  asymptotes  to  its  Tollmi(m-Schlichting  wave  value  at  sufficiently  large  values 
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of  Q* .  We  conclude  that  wall  heating  has  a  progressively  increasing  destabilizing  effect  on 
Tolhnien-Schlichting  waves. 

If  the  wall  is  instead  cooled,  so  that  the  Grasliof  number  is  negative,  then  the  stable 
band  of  modes  in  the  interval  0  <  fl*  <  2.298  for  Go  =  0  increases  with  the  neutral  value 
of  /}*  given  by  (3*\Go\^  when  Go  — >  —  oc.  However  we  again  note  that  for  any  given  value 
of  Go  the  growth  rate  approaches  the  Tolhnicn-Schlichting  value  at  large  values  of  12*. 

8  Conclusion 

We  shall  firstly  consider  the  conclusions  to  be  drawn  about  the  generation  of  vortex 
structures  by  surface  imperfections  when  the  spanwise  lengthscale  is  comparable  to,  but 
shorter  than,  the  body  lengthscale.  We  have  demonstrated  how  the  vortices  develop  in 
a  nonparallel  manner  and  shown  that  a  uniciue  growth  rate  does  not  exist  for  a  growing 
vortex  structure.  We  have  also  shown  the  neutral  curve  associated  with  a  particular  flow 
property  depends  on  the  upstream  history  of  the  disturbance. 

In  Section  4  we  have  discussed  the  localized  forcing  problem  when  the  forcing  operates 
on  a  short  stremwise  lengthscale.  It  was  shown  for  a  forcing  function  of  the  type  considered 
in  Section  4  that  the  forcing  has  the  effect  of  producing  a  similarity  solution  of  the  linear 
disturbance  equations  in  the  region  where  the  forcing  is  applied.  The  similarity  solution 
can  then  be  used  to  form  a  composite  disturbance  field  associated  with  an  isolated  forcing 
function,  if  the  Grashof  number  is  then  varied  we  can  determine  the  effect  of  the  location 
of  isolated  forcing  on  the  onset  of  instability.  Our  results  show  that  there  is  an  optimum 
position  for  the  forcing  which  will  produce  instability  at  the  lowest  value  of  the  local 
Grashof  number,  this  lowest  value  is  about  2. 
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In  Section  5  we  have  considered  the  freestreain  receptivity  problem.  We  have  demon¬ 
strated  the  growth  of  vortex  structures  downstream  and  have  shown  that  the  receptivity 
calculatioi:i  with  b  =  0  leads  to  the  most  dangerous  mode.  In  this  case  Uc  ~  cosaz  at  the 
leading  edge  of  the  wall.  The  explanation  of  this  may  be  that  if  the  disturbance  develops 
in  some  type  of  quasi-parallel  fashion  then  for  higher  values  of  b  the  incoming  distur¬ 
bances  stimulate  the  higher  modes  and  which  arc  more  stable.  It  is  of  interest  to  note 
that  freestream  disturbances  provoke  instability  at  a  much  lower  Grashof  number  than  do 
roughness  induced  motions.  It  would  appear  then  that,  in  an  experiment  where  care  has 
been  taken  to  reduce  the  size  of  disturbances  from  all  sources,  it  will  be  the  freestream 
ones  which  cause  the  growth  of  streamwise  vortices. 

For  the  case  where  the  forcing  varies  on  the  body  lengthscale  as  discussed  in  Section  G 
we  have  demonstrated  how  the  initial  disturbance  decays  and  is  formed  into  a  wake  before 
its  subsequent  rcampliflcation  within  a  wedge  shape  region  further  downstream.  We  have 
noted  that  this  type  of  flow  structure  has  ]:)een  observed  experimentally  for  the  related 
Gortler  type  vortex  problem  of  flow  over  a  concave  wall.  In  fact  the  results  of  Gilpin, 
Imura  and  Cheng(1978)  are  also  consistent  with  this  picture. 

The  results  found  in  Section  7  show  an  unexpected  coupling  between  Tollmicn-Schlichingl 
waves  and  streamwise  vortices  at  low  spanwise  wavenumbers;  in  fact  in  that  regime  the 
two  types  of  disturbances  are  virtually  indistinguishable.  Though  disturbances  with  much 
higher  growth  rates  are  possiljle  at  high  Grashof  numbers,  these  low  wavenumber  distur¬ 
bances  might  be  particularly  relevant  when  the  forcing  mechanism  which  generates  the 
vortices  operates  on  a  long  spanwise  scale.  Moreover  it  could  well  be  that,  even  though 
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larger  linear  growth  rates  are  possible  downstream  where  the  local  wavenumber  has  be¬ 


come  0(Gi4),  the  disturbances  might  be  sufficiently  amplified  near  the  left  hand  branch 
of  the  neutral  curve  for  nonlinear  effects  to  come  into  play.  In  that  case  the  fact  that  larger 
linear  growth  rates  were  available  downstream  would  be  irrelevant. 

We  now  make  some  further  comparisons  with  previous  theoretical  and  experimental 
results.  In  order  to  make  such  comparisons  it  is  convenient  to  define  the  parameters: 

,  Re,  = 

-3 

The  local  Gx'ashof  number  in  our  notation, is  then  given  by  Gx  =  GvxRcx  ^  .  Thus 
if  the  instability  is  caused  by  wall  forcing  we  expect  that  instability  will  occur  whenever 


Gr,Re,~  >~  2.  Wu  an  Cheng  (1976)  made  a  parallel  flow  stability  analysis  of  the 
problem  investigated  here  and  found  that  for  air  instability  occurred  for  GrxRe,~  >~ 
292.  Later  results  given  by  Moutsoglou,  Chen  and  Cheng  (IDSl)  contradicted  those  of  Wu 
and  Cheng  and  Figure  1  of  their  paper  suggest  instability  at  zero  Grashof  number.  This 
result  is  not  unlike  some  of  the  physically  unrealistic  results  given  by  parallel  flow  theories 
of  Gdrtler  vortex  growth.  In  order  to  remove  this  difficulty  Moutsoglou,  Chen  and  Cheng 
retained  higher  order  buoyancy  effects  even  though  they  are  formally  negligible.  The  latter 
approach  is  equivlent  to  the  attempts  made  to  alleviate  the  corresponding  Gdrtler  problem 
by  retaining  higher  order  curvature  effects.  Our  results  show  that  if  nonparallel  effects  are 
accounted  for  in  a  self-consistent  matter  then  instability  occurs  at  a  finite  Grashof  number 
and  the  difficulty  is  not  present. 

Wang  (1982)  investigated  experimentally  the  onset  of  the  vortex  instability,  his  results 

-5 


suggest  that  instability  occurs  for  Gr^Re,  »  >~  55..  This  is  not  consistent  with  our 
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prediction  which  has  the  |  j^ower  replaced  by  |.  We  presume  that  the  experimental  result 
is  not  consistent  with  ours  because  of  the  uncertainty  associated  with  identifying  the  onset 
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Figure  (3.3a-d)  The  downstream  development  of  u,v,w,6  for  the  initial  disturbance  given 
by  (3.7)  with  G  =  .025,  a  =  .069. 


Rg(3.4) 


Figure  (3.4)  The  neutral  curves  for  different  values  of  x  for  the  initial  conditions  (3.7). 


ng(4.3a) 


Fig(4.3b) 


n 

Figure  (4.3a-c)  The  downst.reain  development,  of  tlie  hmrtions  ri,r\w 
(4.22),  the  curves  shown  correspond  to  x=:. 505, .905, 8. 
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Fig(4Jc) 


n 

Figure  (4.3a-c)  The  downstream  development  of  the  functions  u,v,w  associated  with 
(4.22),  the  curves  shown  correspond  to  x=. 505, .905, 8. 


Fig(4.4) 


Figure  (4.4)  The  neutral  curves  corresponding  to  the  downstream  development  shown  in 
Figure  (4.3). 
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(4.23),  the  curves  shown  correspond  to  x=. 505, .905, 8. 
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Fig(4.5c) 


»> 

Figure  (4.5a-c)  The  downstream  development  of  the  functions  u,v,u)  associated  with 
(4.23),  the  curves  shown  correspond  to  x=. 505, .905, 8. 


Fig(4.6) 


Figure  (4.6)  The  neutral  curves  corresponding  to  the  downstream  development  shown  in 
Figure  (4.5). 
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Figure  (5.1)  The  downstream  development  of  u,v,9  for  the  freestream  receptivity  problem. 
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